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5_i ' Abstract 

\ For theoretically consistent determination of a s from jet rates in deep inelastic scattering the 

' dependence on a s of parton distribution functions is in principle as important as that of hard 

scattering cross-sections. For the kinematical region accessible at HERA we investigate in detail 
numerical importance of these two sources of the a s dependence of jet rates. 

> : 

O ■ 1 Introduction 

m ; 

As QCD enters the stage of a truly quantitative theory, great emphasize is put on precise deter- 
mination of its basic parameter, a s , or equivalently A^ , using different sorts of hard scattering 
<3\ ' processes in broad range of kinematical conditions. As a result of numerous phenomenological 

analyses crucial property of QCD - the running of a s - has been verified beyond doubt (for review 
sec ]l|j2[]). However, the increasing precision of these measurements has lead to a discrepancy be- 
tween a s (Mz) extrapolated from measurements of a s (/i) at low scales /j,, and a s [Mz) measured 
directly at LEP. This effect, if confirmed by more accurate data, would indicate that the running 
of a s slows down in the LEP range. This, in turn, might signal the long sought new physics |J. 
. ^ | One of the problems of quantitative determination of the running of a s is related to the fact that 

^ ■ it usually requires combining results of different experiments in different kinematical regions. One 

of the hard scattering processes where the running of ct s (n) can be investigated for the broad range 
of scales in a single experiment is the deep inelastic scattering (DIS) at HERA. The theoretically 
cleanest way for such an investigation, the measurement of scaling violations of F^^x, Q 2 ), is 
marred by the fact that scaling violations are tiny effect at large Q 2 . Recently, however, both the 
HI ||] and ZEUS @ Collaborations have reported evidence for the running of a s (fi) obtained from 
the measurement of jet rates in DIS via the quantity 

p {n 2\ _ ^2+l(Q 2 ) m 

a 1+1 (Q z ) + a 2+ i{Q z ) 

where at+i denotes the cross-section for the production of k hard and one proton remnant jets. 
There are two Monte-Carlo implementations H0 of the NLO calculations of the cross- 



sections ak+i(Q 2 ), which include also the possibility of imposing cuts on various kinematical vari- 
ables. In this paper we shall concentrate on the HI analysis [Q] of 1993 data, which, using PROJET 
4.1 H with the JADE jet algorithm and y c = 0.02, obtained the following result for a s (M^,MS) 

a s (M z ,MS) = 0.123 ± 0.012(stat.) ± 0.008(syst.). (2) 
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This value has been extracted from the Q 2 -dependence of i?2+i(Q 2 ), displayed in Fig. 5 of @, 
using only the two highest Q 2 data points. Within the error bars the value (||) is consistent with 
the world average, but its accuracy is insufficient to draw any firmer conclusions. Discarding (for 
reasons discussed in 0) the three lowest Q 2 points means not taking into account the region 
where most of the running of a s (fi) occurs and where also the data are most sensitive to A^jg . 
The experimental errors would have to decrease by at least a factor of three to make this way of 
determining a s (M^,MS) competitive with the most precise determinations of a s (Mz, MS) based 
on DIS scattering or LEP event shapes data. Some improvement may be expected already from 
the ongoing analyses of 1994 data. 

There is one aspect of the extraction of a s from the measured ratio ([l]) crucial for its consistency. 
In j| a s (fi, MS) (or, equivalently, Aj^g) was considered as a free parameter in the hard scattering 
cross-sections, but not in the parton distribution functions (PDF), for which the MRSH set was 
used. The authors of concede that "for strict consistency" one should "reevolve the distributions" 
using the value of a s (/U, MS) obtained in their paper, but claim that given the presently obtainable 
accuracy "this is not necessary". However, reevolving PDF, that is solving the evolution equations 
with the initial conditions corresponding to the MRSH (or any other chosen set of PDF), but 
using the value of Ay^ extracted from ([[]) instead of the "input" A^jg(MRSH) would solve nothing. 
Except for the case that the extracted Aj^g coincided with the input one, the reevolved PDF 
would be different from the original MRSH ones. This would require iterating the whole procedure 
until the extracted value of A^g coincided with the input one. But then one would have to check 
whether the resulting new PDF still described the whole body of hard scattering data on which the 
original global fit was based. Moreover, although mathematically legal exercise, reevolving PDF 
from unchanged initial conditions with changed A^jg gued in Section 3 below, theoretically 

inconsistent. In the situation when the measurement of the quantity (||) alone is clearly insufficient 
to determine a s , the only meaningful way of using its results is to include them in the conventional 
global fits and see what comes out. Q 

Nevertheless, it might be useful to know under which conditions neglecting A^jg in the PDF 
is numerically good approximation for determining a s from (|l|). The purpose of our paper is to 
investigate this question in detail. The paper is organized as follows. In the next section we recall 
basic relevant formulae and emphasize the importance of considering Aj^g in PDF. In Section 3 a 
simple procedure of constructing PDF for arbitrary Aj^g is outlined. Results of numerical simula- 
tions in which A^jg is varied separately in hard scattering cross-sections and PDF are presented 
and discussed in Section 4, followed by conclusions in Section 5. Detailed elaboration of arguments 
supporting the claim made in Section 2 is left for the Appendix. 



2 Basic formulae 

We use the couplant a = a s (^)/7r as QCD expansion parameter. It satisfies the equation 
da(/i, RS) 



dln^x 



-ba 2 {fi, RS) (l + ca{fi, RS) + c 2 a 2 (/i, RS) + • • •) . (3) 



For 3 colours and rif massless quark flavours b = (33 — 2nj)/6, c = (153 — 19ny)/(66 — 4nj), 
while all higher order coefficients Ck ; k > 2 in (|3|) are free parameters, which define the so called 
renormalization convention (RC). Together with the specification of the initial condition on the 
solution of (H|) they define the renormalization scheme (RS). To get a unique result of perturbative 

2 This is what the MRS group did in their latest global fit |i~2| , which includes beyond the standard set of hard 
scattering data also the very recent CDF results on jet production Ii3[ . 
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calculations both the scale \i and the RS must be specified. Throughout this paper we stay in 
the conventional MS RS and therefore the corresponding specification "MS" in A^jg as well as in 
a s (iJ,, MS) will be dropped in the following. The jet cross-sections at+i in (0) ar e given as integrals 
over the parton level cross-sections Ck+n 

a k+1 {Q 2 ,y c ,A) = V f dxfi(x, M, A)C k+ i t i(Q, M, x, y c , A), (4) 
I Jo 

where fi(x,M) are PDF of the proton evaluated at the factorization scale M and the sum runs 
over all parton species, i.e. quark, antiquaries and gluon. In (||) we have written out explicitly the 
dependence of both jet and parton level hard scattering cross-sections cr k +i, C k+ \^ as well as of the 
PDF fi on A . This dependence reflects the fact that physical quantities do depend on the value of 
A (in a chosen RS). Contrary to that, jet cross-sections u k+ \ evaluated to all orders of perturbation 
theory are independent of both the factorization scale M and the RS. In perturbative QCD the 
parton level hard scattering cross-sections are given as expansions in the couplant a (///A), taken 
at the hard scattering scattering scale f] \i 

C2+i,i(Q,M,x,y c , A) = a( l x/A)\c^l 1>i (x,y c ) + a( t x/A)cfl lji (Q/M, x,y c ) \ , (5) 
C 1+1>i (Q,M,x,y c ,A) = <£l u (Q,x) + a(ji/A)<j& u (Q/M,x,y e ). (6) 



In standard global analyses of hard scattering processes ]B|-|16|] , A is fitted together with a set of 
parameters af , describing distribution functions p^\x, M) of parton species j at some initial scale 
Mq, usually in the form 



p®(x, M ) = afx^ (1 - xf^ (l + af^ + afx + afx^ 2 ) . (7) 



To determine the relative importance of varying A in PDF and hard scattering cross-sections, we 
write the derivative da"2+i/dlnA: 

da 2+1 (Q 2 ,A) _^ fi A _ r d/<(s,M,A) a(M) ^ +/i(g>M>A)<g (g) da(M) 



dlnA 



dlnA r ' w ' ,r ' ..v-. v-/ dlnA 



a ( M )E [ dx [-[ jfj(y,M,A)P^(z)a(M)4l hi (x) + bf i (x,M,A)a(^^ 



,(8) 



where z = x/y and P-p(z) are the LO AP branching functions []. In the above expression we have 
used @ and standard evolution equations for PDF, exploited the equality da(M/A)/dlnM = 
— da(M/A)/dln A and took into account that the same relation between derivatives with respect 
to M and A holds also for PDF. The first term in the brackets of (||) results from variation of A in 
PDF while the second from variation of A in parton level hard scattering cross-sections. Eq. (||) 
implies that the leading order term of do"2+i/dlnA involves merely the LO evolution equations for 
PDF and leading order term of parton level hard scattering cross-sections C^+n- The dependence 
of (T2+1 ° n A starts therefore already at the LO and consequently also the relative importance of 
the two terms in the brackets of (|8|) is basically the LO effect. Note, however, that associating 
A in o"2+i with any well-defined RS requires including the NLO corrections as well. Eq. (||) tells 
us that the variation of A in the couplant a(/i/A) changes the jet cross-section 02+1 by the term 

''Although in general fj, 7^ M, we shall in the rest of this paper follow the usual practice and set \i = M. 
4 In (8) we have for brevity suppressed the dependence of c^l-y ; on y c and written a(M/A) simply as a(M). 
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of the same order (schematically 0{a)a2+\) as when A is varied in PDF. Consequently if we want 
to determine A from jet cross-section 02+1, we must vary it in PDF as well. 

For the quantity i?2+i(Q 2 ) the situation is less obvious. PDF appear in both the numerator and 
denominator of (|l|) and so some cancellations might occur, while the hard scattering cross-sections 
start as 0{a) for C2+1 i and as 0(1) for C14.1 j. Nevertheless as shown in the Appendix also for the 
ratio R2+1 varying A in PDF generates terms which are of the same order as those resulting from 
the variation of A in the hard scattering cross-sections Ck+n- 



3 Parton distribution functions for arbitrary A 

In the previous Section we have shown that for a theoretically consistent extraction of A from 
i?2+i(Q 2 ), A mus t be varied not only in hard scattering cross-sections, but also in the PDF. We 
now want to investigate the numerical importance of varying A in PDF. For this we first have 
to know what to do with the parameters Mo,aj , specifying the initial conditions @, when A is 
varied. The obvious choice seems to be to keep them fixed and simply solve the evolution equations 
with variable A . This, however, would be inconsistent with the idea of factorization, which implies 
that also the initial conditions depend on A! 

To see the caveat, let us first consider as well as the evolution equation for the nonsinglet 
quark distribution function q^s(x,M) at the LO. In term of moments, defined as 



/(«) 



we have 



wherefrom 



dg m (n,M, A) 
dlnM 



dxx n f(x) 



a(M/A)q m (n, M)P® (n) , d n = -P^> (n)/b 



(0), 



<ZNs(n, M, A) = An [a(M/A)} d " , o(M/A) 



In the above expressions Pj^g (n, M) are moments of LO nonsinglet AP branching function and 
A n are dimensionless constants describing the nonperturbative, long distance, properties of the 
nucleon. These constants are unique pure numbers. Via (11) they determine the asymptotic 
behavior of (/Ns( n > M, A) for M — > 00 JlTj and thus provide alternative way of specifying the initial 
conditions on the solutions of (]To[) . The separation in (|ll]) of q^s(n, M, A) into two parts - one 
calculable in perturbation theory and the other incorporating all the nonperturbative effects - is the 



61n(M/A)' 



0) 



(10) 



(11) 



essence of the factorization idea [17]. Note that in ( |TT| ) the dependence on A is simply a reflection 
of its dependence on M. Knowing the latter, we know the former. For finite initial Mq ( |TT[) implies 



q m (n,M,A) = g NS (n,M ,A) 



a(M/A) 
a(M /A) 



where 



/ln(M/An 
V ln(M /A)y 



In 



= qns{n, M , A) exp(-d n s) 
o(M /A)' 



(12) 



(13) 



a{M/A) 

is the so called "evolution distance" fig] . The second equality in ( |i~3| ) holds exactly at the leading 
order and approximately at higher ones. Rewriting (^) again in the form (|TTj) 



q^ S (n,M, A) 



g N s(ra,M , A) 
(a(M /A)) d " 



(a(M/A)Y 



(14) 
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we see that the ratio of gNs( n ; Mq, A) and (a(Mo/A)) rfn equals the invariant A n and thus is both 
Mq and A independent. This, however, is possible only if the initial moments q^s(n, Mq, A) do, as 
already indicated, depend on A as well! What happens in (14) is that the A and Mq dependen- 
cies coming from a(Mo/A) in the brackets of cancel those of the initial q^s( n -, A, Mq). Note 
that considering q^s( n , M, A) as a function of s, formula ( |i~4| ) implies that entire dependence of 
5Ns( n ) M, A) on A comes solely from the term ln(M/A) = l/a(M/A) in (13)! The same happens 
in the realistic case of coupled quark and gluon evolution equations |l8| . The dependence of PDF 
on M and A through the evolution distance variable s is then a direct consequence of the above 
relations for the moments. 

To investigate the dependence of PDF on A we therefore cannot fix initial conditions at some 
finite Mq, specified by p^\x, Mq) in (Q), and vary A in the evolution equations alone. The preceding 
discussion also suggests what can be held fixed in the process of varying A - the constants A n \ 
There is no principal obstacle to carrying out QCD analyses with the initial conditions specified 
by means of the invariants A n , which together with A should then be fitted to data. There is also 
a well developed technique of orthogonal polynomials expansions fl^| , which could be used for this 
purpose. There are, however, technical problems in practical implementation of this method. One 
of them concerns the number of parameters A n required for a specification of initial conditions. 
As shown in pOfl , one needs at least 9 to 10 first moments of PDF to achieve sufficient precision 
in reconstruction of PDF from their moments. That is about twice the number of free parameters 
used in conventional initial conditions (0) at finite Mq and seems too much for any reasonable QCD 
fit. 

However, the preceding paragraph also suggests a simple procedure, which makes use of some of 
the available parameterizations and generalizes them to arbitrary A , keeping, as discussed above, 
the constants A n fixed. From all the available parameterizations of PDF those given by analytic 
expressions of the coefficients a[ on a general factorization scale M via the variable s, defined 
in (|l3|), are particularly suitable for this purpose. In our studies we took several of the GRV 
parameterizations, both LO and NLO ones, defined in At the LO the recipe for construction 
of PDF for arbitrary A is simple: 

• Use any of the parameterizations expressing the solutions of evolution equations as functions 
of the "evolution distance" variable s, defined in (|l3|). 

• Keep A in ln(Mo/A) fixed at the value obtained in the fit. 

• Vary A in ln(M/A). 

• Use the original parameterization with variable s = s(A). 

This construction is exact at the LO. At the NLO the couplant a(M/A) is no longer given simply as 
l/ln(M/A) and consequently the second equality in ( |l3| ) holds only approximately. This, together 
with the additional term appearing in the NLO expression for (?Ns( n ) M, A) 



q NS (n,M,A,A n ) = A n 



a(M/A) 



-i d n 



1 + ca(M/A) 



l + ca(M/A)) P Ns (n)/6c , (15) 



destroy simple dependence of </ns on M, Mq and A via the variable s. Using the above recipe for 
the NLO GRV parameterizations provides therefore merely an approximate solution of our task, 
but one that is sufficient for the purposes of determining the relative importance of varying A in 
PDF and hard scattering cross-sections. 

As an illustration of our procedure we plot in Fig. 1, for x > 10~ 4 and two values of Q 2 , 
the dependence of valence and sea it-quark and gluon distribution functions on A. The curves 
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correspond to the LO GRV parameterization of [21]. We see that this dependence decreases with 
increasing Q 2 , is most pronounced for the gluon distribution function and almost irrelevant for the 
valence quark one. These features are qualitatively the same as those obtained in refs. [22j,23|, which 
contain results of global fits performed for several fixed values of A. Nevertheless, quantitatively 
the sensitivity to the variation of A , observed in these papers, is markedly weaker than that of our 
results, displayed in Fig. 1. The reason for this difference is simple. In global analyses of [22, 2^] 



fitting the data for different values of A leads to different values of boundary condition parameters 
, which partially compensates the variation of A. Eqs. (pd|) or (15) tell us that for any fixed 



scale M and for any moment n (and thus for any x) the change of A can be fully compensated by 
an appropriate change of the constants A n . For a finite interval of M this compensation cannot be 
complete, but it is clear that the dependence of the fitted PDF on A comes out weaker than in our 
procedure where only A is varied and all the coefficients A n (or, analogously, af ) are held fixed. 

4 Numerical results 

All the results reported below were obtained using the PRO JET 4.1 generator ||. To quantify the 
dependence of -R2+1 on A separately in hard scattering cross-sections and PDF we consider it as a 
function of Q 2 ,y c and two independent A-parameters, denoted as A hsc and A pdf respectively. The 
results are then studied subsequently as a function of A hsc , A pdf , Q 2 and y c . First we analyse the 
case with no cuts on the final state jets and then discuss the changes caused by the imposition of 
cuts as specified in Q. 

4.1 The case of no cuts 

The simulations were done at both the leading and next-to-leading order and for 

• 12 values of Q 2 (equidistant in lnQ 2 ) between 5 and 10 4 GeV 2 , 

• 5 values of AP df or A hsc , equal to 0.1,0.2,0.3,0.4 and 0.5 GeV, 

• 3 values of y c = 0.01, 0.02, 0.04. 

First we discuss the results of the NLO analysis. Fig. 2 shows a series of plots corresponding to 
y c = 0.01 and displaying the dependence of R2+1 on A hsc for fixed A pdf (dotted curves) or on 
A pdf for fixed A hsc (solid curves). Each plot corresponds to one of the 6 selected values of Q 2 . 
We see that for Q 2 below about 40 GeV 2 the solid curves are steeper that the dotted ones, while 
above 40 GeV 2 the situation is reversed. This means that below 40 GeV 2 varying A in PDF is 
more important than varying it in the hard scattering cross-sections for Q 2 , while the opposite 
holds above 40 GeV 2 . The same information as in Fig. 2 is displayed in a different manner in 
Fig. 3, where the Q 2 dependence of -R2+1 is plotted for various combinations of A pdf and A hsc . 
The shape of the curves in Fig. 3 is a result of two opposite effects: as Q 2 increases the phase 
space available for produced jets increases as well, but at the same time the value of the running 
a s (Q/A) decreases. At low Q 2 , on the other hand, phase space shrinks, but a s (Q/A) grows. The 
same features of the relative importance of varying A pdf and A hsc as mentioned above are reflected 
in a bigger spread of the five curves in Figs. 3d)-f) compared to those in Fig. 3a)-c) for Q 2 < 40 
GeV 2 and smaller spread above 40 GeV 2 . Note also that for fixed A pdf the sensitivity to the 
variation of A hsc is about the same at all Q 2 . The shape of curves in Fig. 2 suggests that the ratio 

V{Q ,y c ,r,s) = — ^ — m ^ r, (16) 

R 2+ i {Q z , y e , A pdt = r, A nsc = s) 
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is to a good approximation a linear function of s for any fixed Q 2 ,y c ,r. Considered as a function 
of s for fixed r, ( |i~6| ) quantifies the relative importance of varying A in the PDF (described by R2+1 
in the numerator) and in hard scattering cross-sections (described by R2+1 in the denominator). 
To summarize the results of Fig. 2 we fitted V(Q 2 ,y c , r, s) by a linear function of s and defined 

^^ gS^ (17) 

which, by construction, is s-independent function of Q 2 ,y c and r. Positive W means that the 
variation of A in the PDF is more important than that in the hard scattering cross-sections, while 
for negative W the situation is opposite. The Q 2 -dependence of W(Q 2 ,y c , r) is plotted in Fig. 4a-c 
for three values of y c = 0.01, 0.02, 0.04 and five values of r = 0.1, 0.2, 0.3, 0.4, 0.5 GeV. In Fig. 5a we 
compare the Q 2 dependence of W(Q 2 ,y c ,r = 0.2 GeV) for different values of y c = 0.01,0.02,0.04. 
Similar plots can be drawn for other values of r as well. From Figs. 2-5a we conclude that at the 
NLO: 

• For Q 2 below about 40 GeV 2 the variation of A is more important in PDF than in the hard 
scattering cross-sections, while above 40 GeV 2 the situation is reversed. 

• Above Q 2 ~ 10 3 GeV 2 the variation of A in PDF becomes negligible. 

• The preceding conclusions depend only weakly on y c . 

At the LO the general features are the same as those displayed in Figs. 2-5a and we therefore 
merely summarize them in Fig. 5b. Comparing Figs. 5a and 5b we see that at the LO 

• the "cross-over" point Q 2 r , where W(Q 2 r ,y c ,r) = 0, is at higher values of Q 2 r ~ 50 — 200 
GeV 2 and depends more sensitively on y c , than at the NLO, 

• the relative importance of varying A in PDF is for almost all values of Q 2 and y c larger than 
at the NLO, 

• the difference between the LO and NLO results increases with increasing y c . 
4.2 The case of HI cuts 

In [|j two subsamples of events with jets are defined: 

a) high Q 2 sample: Q 2 > 100 GeV 2 , 

b) low Q 2 sample: 10 < Q 2 < 100 GeV 2 . 

In both subsamples outgoing jets angle in the HERA lab system was restricted to lie in the interval 
10° < 0j£ < 145°, but there are also cuts on other variables, in which the two samples differ. In 
Fig. 6 plots analogous to those of Fig. 3, but incorporating the HI cuts, are displayed. We see 
that the experimental cuts suppress the sensitivity of i?2+i(Q 2 ) to the variation of A in the PDF 
in the low Q 2 region, while leaving it essentially unchanged the high Q 2 one. We conclude that in 
the region Q 2 > 100 GeV 2 , used in Q for the extraction of A from the measured i?2+i(Q 2 ), the 
variation of A in the PDF can be safely neglected. 
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5 Conclusions 



In this paper we have shown that theoretically consistent determination of a s from R2+1 in DIS 
requires the variation of A in both the hard scattering cross-sections and parton distribution func- 
tions. After defining a procedure for the construction of PDF which satisfy the corresponding 
evolution equations for an arbitrary value of A, we used them to investigate the numerical impor- 
tance of varying A in PDF for the quantity it^+i- The results of this investigation show that if no 
cuts are applied on produced jets the variation of A in PDF 

• is more important at the LO than at the NLO, 

• is more important than that in the hard scattering cross-sections in low Q 2 region (Q 2 < 40 
GeV 2 at the NLO and Q 2 < 50 - 200 GeV 2 at the LO), 

• can be safely neglected above Q 2 1000 GeV 2 . 

We have also demonstrated that in the low Q 2 region the sensitivity of i?2+i(Q 2 ) to the variation 
of A in PDF is strongly suppressed by imposing experimental cuts applied in the HI paper j|. 
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Appendix 

To show that A in PDF must be varied also in the ratio (Q) we shall draw on analogy between the 
cross-sections cifc+i and the longitudinal and transverse structure functions Fl, F 2 given as 



F 2 (x,Q 2 ,A) fidy 



x y 



12 



12 



£ e 2 qi (y, M, A)Cf > (z, Q, M, A) + £ e 2 G(y, M, A)^ (z, Q, M, A) 



F L (x,Q 2 ,A) = /•• (by r 
x Jx y 



1=1 
12 



]T e 2 qi (y, M, A)Cf > (z, Q, M, A) + £ e 2 G(y, M, A)C ( G L) (z, Q, M, A) 

_i=l i=l 

where z = x/y, the sums run over all quark and antiquark flavours and 

C®(z,Q,M,A) = 5{l-z) + a{M/A)cf\z,Q/M)+0{a 2 {M/A)), 

C { c\z,Q,M,A) = a{M/A)c§(z,Q/M) + 0{a 2 {M/A)), 

C( q L \z,Q,M,A) = a{M/K)cf\z) + 0{a 2 {M/K)) 1 

C^ ] {z,Q,M,A) = a{M/A)c^\z) + 0{a 2 {M/A)). 

The quantity analogous to R2+1 is the ratio 

F L (x,Q 2 ,A) 



R L (x,Q 2 ,A) = 



In terms of moments (p~8[)-(19) read 



12 



F 2 (x,Q 2 ,A)+F L (x,Q 2 ,A)' 



12 



(18) 

J 

(19) 

(20) 
(21) 
(22) 
(23) 

(24) 



F 2 (n - 1, Q 2 ,A) = $*(n, M, A)C< 2 >(n, Q, M, A) + ^ e 2 G(n, M, A)cg\n, Q, M, A), (25) 



i=l 
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12 12 

F L (n - 1, Q 2 ,A) = £ e 2 iqi (n, M, A)Cf >(n, Q, M, A) + £ e 2 G(n, M, A)4 L) (n, Q, M, A). (26) 
i=i i=i 

For moments the quantity analogous to i?2+i(Q 2 ) is the ratio 

75, ^2 ^_ F L (n-l,Q 2 ,A) 

i? L (n,Q ,A = — TT^TT- ( 27 ) 

F 2 (n - 1,Q 2 ,A) + F L (n - 1,Q 2 , A) 



Note that i?i(n, Q 2 , A) defined in ( |2"7| ) is not the moment of the ratio Rl(x), but rather the ratio of 
the moments! We mention the analogy of i?2+i to Rl(x) and i?i(n) because all three quantities have 
similar structure, but the later two are simpler and, moreover, Rl{u) can be treated analytically. 
The technique of moments has recently been applied in [24 1 to a related problem of extracting gluon 
distribution function from jet production in DIS. 

In the expressions for ak+i, R2+1, as well as for F2 ) l,Rl{x) an d Rl, we have written out 
explicitly their dependence on the value of A. Contrary to that, these observables are formally (i.e. 
order by order of perturbation theory) independent of the factorization scale M, the M-dependence 
of PDF being cancelled by that of the hard scattering cross-sections. In the case of Rl(u,Q 2 , A), 
defined in (|27|), the LO expression reads 



R L (n,Q*,A) = -(M/A) "Iff' W + (S C(". W 



a(M/A) 



(E,ef)G(n,M,A)cg ) (n) 
q E t eMn,M,A) 



+ •■-. (28) 



We see that in the quark contribution to the LO term in ( p7| ) quark distribution functions in the 
numerator and denominator of (^) cancel out 

i?r rkS = «(M/A)4 i )+... (29) 

and R L is thus equal to the parton level hard scattering cross-section. This, however, is 
not the case for the second term in d28|), corresponding to gluon contribution, which behaves 
rather like a l+ ^ n \ where A(n) is, roughly speaking, the difference between gluon and quark LO 
anomalous dimensions f| The cancellation of quark distribution functions in the LO coefficient 
of (27) is a consequence of the fact that convolutions appearing in (|i~8|) and (|l9|) translate to 



simple multiplications for the respective moments. The above mentioned cancellation holds only 
for the ratio (27) of moments but not for the ratio ( |24| ) of structure functions themselves. This is 



obvious when we write i^2,i(a?, Q 2 ) as expansions in orthogonal polynomials, which contain infinite 
sums over different moments of PDF p9| , [20[| . As a result, even in the quark contribution to -R2+1 
variation of A in PDF is in principle equally important as that in hard scattering cross-sections. 
This conclusion holds for the ratio i?2+i(Q 2 > A) as well. 
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Figure 1: The dependence of xu v (x, Q 2 ), xu sca (x, Q 2 ) and xG(x,Q 2 ) of the value A pdf according 
to the recipe described in the text and for two values of Q 2 . At low x, the curves correspond from 
above to the values of A pdf , subsequently, 0.1,0.2,0.3,0.4 and 0.5. 
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Figure 2: The ratio R 2 +i(Q, y c , A pdf ,A hsc ) as a function of A pdf for fixed A hsc and a series 
of Q 2 values (solid lines), or as a function of A hsc for fixed A pdf and the same set of Q 2 values 
(dotted lines). The dashed lines correspond to simultaneous variation of A pdf = A hsc . Solid curves 
correspond, from below, to fixed A hsc = 0.1, 0.2, 0.3, 0.4, 0.5 and the dotted to the same fixed values 
of A pdf . In all plots y c equals 0.01. 
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Figure 3: (a-c) Q 2 -dependence of the ratio R 2 +i(Q, y c , A pdf , A hsc ) for fixed AP df = 0.1,0.3,0.5 
GeV and five values of A hsc = 0.1,0.2,0.3,0.4,0,5 GeV. (d-f) The role of A pdf and A hsc is 
reversed. No cuts were applied in evaluating R2+1 and y c = 0.01. In all plots the curves are 
ordered from below according to increasing A hsc (in a-c) or A pdf (in d-f). 
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Figure 4: (a-c) The quantity W(Q 2 ,y c ,r) as a function of Q 2 for five values of r (given in GeV) 
and y c = 0.01,0.02,0.04. All curves correspond to the case of no cuts. 
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Figure 5: The Q 2 -dependence of W(Q 2 ,y c ,r = 0.2) for three different values of y c at the NLO (a) 
and LO (b). All curves correspond to the case of no cuts. 
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Figure 6: The same as in Fig. 3 but for the HI acceptance cuts. 
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